Proc. 20th Winter Workshop on 
Nuclear Dynamics (2003) 000-000 



20th Winter Workshop 
on Nuclear Dynamics 

TVelawny Beach, Jamaica 
March 15-20, 2004 



Hadronization via Coalescence 
V. Greco 1 and CM. Ko 1 



o 
o 

Cyclotron Institute and Physics Department, Texas A&M University, 
TX-77843-3366 College Station, USA 

Abstract. We review the quark coalescence model for hadronization in rel- 
ativistic heavy ion collisions and show how it can explain the observed large 
baryon to meson ratio at intermediate transverse momentum and scaling of 
the elliptic flows of identified hadrons. We also show its predictions on higher- 
■ order anisotropic flows and discuss how quark coalescence applied to open- 

and hidden-charm mesons can give insight to charm quark interactions in the 
ly-j ■ quark-gluon plasma and J/^ production in heavy ion collisions. 
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O \ 1. Introduction 

To find the signals of the quark-gluon plasma and to study its properties in rela- 
tivistic heavy ion collisions, it is important to understand how quarks and gluons 
are converted to hadrons during hadronization. Because of its non-perturbative 
nature, hadronization has so far been treated only phenomenologically based on 
the statistical model, the duality model, or the coalescence model. The coalescence 
model, in which colored partons are combined into singlet hadron clusters, was first 
suggested as a possible mechanism of hadronization in studying pion correlations 
in relativistic heavy ion collisions [ ^ . It was later applied to describe the chem- 
ical composition of the hadronic matter produced in heavy ion collisions at SPS 
and RHIC [ 2 . More recently, the quark coalescence model has been shown to 
explain not only the larger baryon to pion ratio at intermediate transverse momen- 
tum but also the scaling relation among the elliptic flows of identified hadrons that 
are observed in Au+Au collisions at the Relativistic Heavy Ion Collider (RHIC) [ 
01 El El El 13 El- In this talk, the quark coalescence model and its applications to 
heavy ion collisions at RHIC will be reviewed. Besides particle spectra and elliptic 
flows, we also discuss results from the coalescence model on resonance production, 
higher-order V4 flow, and charmed hadron production and flow. 
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2. The coalescence model 



In the coalescence model, the transverse momentum spectrum of hadrons that con- 
sist of n quarks is given by the overlap of the product of n quark phase-space distri- 
bution function f q (xi,Pi) with the Wigner distribution function fn(xi..x n \pi..p n ) 
of the hadron, multiplied by the probability gn of forming from n colored quarks a 
color neutral object with the spin and isospin of the hadron, i.e., 

§r^=9H J J! P '( 2 t)3^f l f q (xi,Pi)fH(xi..x n ; Pl ..p n ) (pt - ' 

where da denotes an element of a space-like hypersurface. 

The n quark phase space distribution is usually approximated by the product 
of the single quark distribution function, which consists of both a thermal and 
a minijet component, separated by the transverse momentum pq ~ 2 GeV. For 
heavy ion collisions at RHIC, the quark transverse momentum spectrum below pq 
is taken to be thermal with a temperature T = 170 MeV and a radial flow velocity 
/3 = 0.5, which are consistent with both experimental data and predictions from 
hydrodynamical calculations. The masses of thermal quarks are taken to be those 
of constituent quarks, i.e., m U;< i = 300 MeV and m s = 475 MeV. Above po, partons 
are from the quenched pQCD minijets with power-law like transverse momentum 
spectrum [[5], and their masses are those of current quarks. Both soft thermal and 
hard minijet partons are assumed to distribute uniformly in a fireball of transverse 
radius of 8.15 fm and longitudinal length of 4.35 fm, corresponding to a volume of 
900 fm 3 . The number of quarks are then fixed by the measured transverse energy 
of about 700 GeV for Au+Au collisions at ^/JWn = 200 GeV. 

For the Wigner distribution functions of hadrons, they are taken to be a sphere 
in both space and momentum with radii A r and A p , respectively, which are further 
related by A r • A" 1 = 1 according to the uncertainty principle. A good description of 
both pion and proton spectra can be obtained with a radius parameter of A p = 0.24 
GeV for mesons and 0.35 GeV for baryons. 



3. Particle spectra and elliptic flows 

In coalescence model, hadrons are formed from quarks that are close in phase space. 
As a result, baryons with momentum P t are produced from quarks with momenta 
~ Pt/3, while mesons with same momentum are from quarks with momenta ~ Pt/2. 
Since the transverse momentum spectra of quarks decrease with px, production of 
high momentum baryons from quark coalescence is more favored than mesons. This 
is opposite to the fragmentation process where baryon production is penalized with 
respect to mesons as more quarks need to be produced from the vacuum, leading 
to a typical p/n ratio of ~ 0.2. Results on p/n ratio based on Eq.JI} are shown 
in FiglUleft) together with data from PHENIX [ ITH]. Although different methods 
have been used in evaluating the coalescence integral in Eq.JQ (for a review see [ 
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111)'), they all lead to enhanced p/ir ratio (and similarly the K/A ratio). As shown in 
Pignjleft), our approach [El EI) which includes resonance decays (to be described 
in the next paragraph) and avoids the collinear approximation by using the Monte 
Carlo method to evaluate the multi-dimensional coalescence integral, also gives a 
good description of the p/ir ratio at pr < 2 GeV. 
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Fig. 1. Left panel: Antiproton to pion ratio in Au+Au collisions at v /snn=200 

GeV. Solid and dashed curves are results with and without contribution to antipro- 
ton from coalescence of thermal with minijet partons. Filled squares are experi- 
mental data from PHENIX [EJ- Right panel: Elliptic flows of identified hadrons. 
Lines are from the coalescence model, while symbols are data from STAR [E] and 
PHENIX [|H|. 

Another typical feature of quark coalescence is the scaling of hadron elliptic 
flows according to the number of constituent quarks in a hadron: V2 J H{PT)/n — 
v 2,q{pr /n), which can be derived under the approximation that only quarks with 
equal momentum can coalescence into hadrons [ 8 . Without any constraint on the 
relative momentum of coalescing quarks, the scaling relation between baryon and 
meson elliptic flows can be badly violated [E]- In our coalescence model, we have 
taken into account the finite momentum distribution of quarks inside hadrons, and 
this leads to only a small breaking (~ 5%) of the so-called "coalescence scaling" [ 
E], which is seen to hold within the errors of experimental data [ 1131 114j . 

In Fig^ (right), we show how the results from our coalescence model can repro- 
duce the available experimental data [El El for tt, K,p, A once the quark elliptic 
flow is extracted from a fit to the pion data [El- We note that the contribution of 
hadrons from minijet fragmentation is not included in this calculation. Its inclusion 
would lead to a universal hadron elliptic flow at momentum above pr ~ 6 GeV [ 
7 . In Ref.[E], the quark elliptic flow was extracted from a fit to the pion data, 
and its shape, especially the saturation at px ~ 1 GeV, is in agreement with that 
from parton cascade calculations [I15j. At low momenta, the elliptic flow from the 
hydrodynamical model at the phase transition temperature could also be used in 
the coalescence model calculations for hadron elliptic flows. This is because the 
description of elliptic flows at low pr are similar in the hydrodynamical and the 
coalescence model, once we correct for energy conservation in quark coalescence 
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and include the residual radial flow effect from the hadronic stage. 

Another possible source of scaling breaking in hadron elliptic flows is the feed- 
down from resonance decays. We have studied how elliptic flow of stable hadrons 
are affected by contributions from the decay of resonances [ E| • Particles from 
these decays have in principle a different elliptic flow from that of resonances, and 
its value depends on the competition between a shift in their momenta, which gives 
a larger elliptic flow, and the randomization of their momenta, that decreases their 
initial azimuthal anisotropy. ft turns out that particles like p, A, and K from 
resonance decays have elliptic flows that are very similar to the directly produced 
ones. Therefore, the inclusion of resonance effect does not destroy the coalescence 
scaling of these stable hadrons. On the other hand, pions from the decay of ui, K*, 
A show a significant enhancement of their elliptic flow at pt < 2 GeV. Although, the 
effect of these resonances is reduced by pions from rho meson decays, which have 
an elliptic flow more close to that of direct pions, the overall effect of resonance 
decays on pion elliptic flow is still non-negligible. The breaking of coalescence 
scaling due to resonance decays together with that due to finite quark momentum 
spread mentioned above lead to a better agreement with available data as shown 
in Fig^left), where the decay of resonances was included [El Ej- However, one 
should be aware that such an agreement may be fortuitous especially at very low 
Pt, again because of the lack of energy conservation in quark coalescence. 

Higher-order momentum anisotropy have also be measured recently in experi- 
ments. In particular, the fourth-order harmonic V4 was found to be sizeable [I17|. as 
first suggested in Ref.[ll8|. Higher-order flows provide the possibility to further test 
the coalescence picture for hadronization. Based on the naive coalescence model 
that only allows quarks of equal momentum to form a hadron, it has been shown 
that the meson V4 m and baryon V4 b are related to quark V2,q and v^ q by [I19|: 

V4,m(2Pt) ^ 2t>4 : q + v\ q v 4 , b (3pt) ^ 3u 4 , 9 + 3^, (2) 

with u n ,g's evaluated at pt- We can see that the ratio between baryon and meson 
Vi is 3 if there is no Vi, q at quark level. However, a fit to the charged particle v 4 
data using the quark coalescence model requires that v^ q ~ 1.8w| q , which is similar 
to «4, g ~ w|,g f rom the parton cascade [E3- With non-zero v^ q , v^b/v^m is thus 
expected to be less than 3. These results are shown in FigO (left) together with 
experimental data on 1*4 of charged hadrons. 

4. J/^ and D meson 

The coalescence model can also be used to study open and hidden charmed hadron 
production in relativistic heavy ion collisions. With the large heavy quark mass, 
a more rigorous treatment of the coalescence process can be achieved [EOl- Study 
of charm meson production at RHIC using the quark coalescence model allows us 
to address the following issues: (a) What is the sensitivity of the J/\P abundance 
and its px spectrum to the momentum distribution of charm quarks? (b) How 
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Fig. 2. (Left) V4 for charged hadrons, pion and proton from quark coalescence 

together with the experimental data (STAR) for charged hadrons (Right) 
Elliptic flow for D mesons and single electrons from their decays. 



does the interplay of charm- and light-quark distributions translate into the ellip- 
tic flow of £)-mesons? Without assuming complete thermalization of charm quark, 
which is common to statistical models, makes it possible to discriminate the rela- 
tive importance between J/tp regeneration and suppression in relativistic heavy ion 
collisions. Studies based on the rate equation have shown that regeneration of J/\l/ 
can bring its abundance back to the initial value expected from the superposition 
of nucleon-nucleon interactions [ I21| . On the other hand, the initial momentum 
distribution of charm quarks from the PYTH1A routines is much broader than if 
they are thermalized in the quark-gluon plasma. The average relative momentum 
between two quarks is thus larger than that of the charm quarks inside the char- 
monium. Quantitatively, this causes a decrease of their recombination probability 
into J/ty by about a factor 3 with respect to the thermal case [I22|. 

Another way of probing charm interactions in the medium is through the elliptic 
flow of single electrons from D meson decays. We have already shown how quark 
coalescence is able to describe the elliptic flow of light hadrons. Our knowledge 
of i>2(pt) for light hadrons thus allows for a quantitative prediction of V2d(pt) for 
the case in which charm quarks do not experience any final-state interactions. This 
represents the lower limit for the V2D of D mesons (D w/o in FigEl (right)), because 
it is only due to the elliptic flow of light quarks. An upper limit can be obtained 
by assuming that charm quarks have a V2D equal to that of light quark (D w in 
Fig'EJright)). It is also seen in Fig^right) that single electrons (diamonds) from 
D meson decays have essentially the same flow as that of D mesons. 

In our calculations, we have not included D mesons produced from jet fragmen- 
tation. The elliptic flow of these D mesons depends on the energy loss of charm 
quarks in the quark-gluon plasma. Since the latter is predicted to be smaller than 
that of light quarks, a quite small i>2 (a first estimation put the upper limit around 
5% [ I23p is expected from these D mesons. Experimental measurement of V2D 
thus allows to distinguish between the quark coalescence and the jet fragmentation 
mechanism for D meson production in relativistic heavy ion collisions. 
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